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The family of Gowdy universes with the spatial topology of a three-torus is studied both classically 
and quantum mechanically. Starting with the Ashtekar formulation of Lorentzian general relativity, 
we introduce a gauge fixing procedure to remove almost all of the non-physical degrees of freedom. 
In this way, we arrive at a reduced model that is subject only to one homogeneous constraint. The 
phase space of this model is described by means of a canonical set of elementary variables. These 
are two real, homogeneous variables and the Fourier coefficients for four real fields that are periodic 
in the angular coordinate which does not correspond to a Killing field of the Gowdy spacetimes. We 
also obtain the explicit expressions for the line element and reduced Hamiltonian. We then proceed 
to quantize the system by representing the elementary variables as linear operators acting on a vector 
space of analytic functionals. The inner product on that space is selected by imposing Lorentzian 
reality conditions. We find the quantum states annihilated by the operator that represents the 
homogeneous constraint of the model and construct with them the Hilbert space of physical states. 
Finally, we derive the general form of the quantum observables of the model. 
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I. INTRODUCTION 

The alternative formalism for general relativity put for- 
ward by Ashtekar [1,2] has renewed the hopes of consis- 
tently quantizing the gravitational interaction in a non- 
perturbative way. In contrast to the situation found in 
the geometrodynamic formulation, the gravitational con- 
straints acquire a simple, polynomic form in terms of the 
Ashtekar canonical variables. Besides, by shifting the em- 
phasis from geometrodynamics to connection dynamics, 
the introduction of the Ashtekar variables has allowed the 
use in gravity of mathematical techniques that had been 
developed in the quantization of gauge field theories. 

In order to gain insight into the kind of problems that 
one will probably have to face when quantizing full gen- 
eral relativity, a lot of attention has been devoted in 
the last years to the quantization of gravitational mod- 
els with different types of spacetime symmetries [3,4]. 
Most of the systems studied are however minisuperspace 
models [3]. These are clearly inadequate to discuss the 
difficulties that will presumably arise in the quantization 
of full gravity owing to the presence of an infinite num- 
ber of degrees of freedom. A possible way to analyze 
such difficulties would be to consider the quantization of 
midisuperspace models. The symmetry of this type of 
models is not large enough as to eliminate all the local 
degrees of freedom, so that their quantization will lead 
to a true quantum field theory. 

In a recent paper [4], Ashtekar and Pierri carried out 
the quantization of the Einstein-Rosen cylindrically sym- 
metric spacetimes [5], completing previous works on the 
subject by Kuchaf [6] and Allen [7]. To our knowledge, 
this is the only gravitational midisuperspace model that 
has been rigorously quantized in the literature. 



It would be of interest to have at our disposal other 
examples of midisuperspace models whose quantization 
can be achieved. Natural candidates for such models are 
provided by spacetimes with two commuting spacelike 
Killing fields [8]. These spacetimes can generally be de- 
scribed by two local physical degrees of freedom which 
depend only on one of the spatial coordinates. Besides, 
since the pioneering work by Geroch [9] , it is known that 
the Einstein equations of these spacetimes present an in- 
finite number of symmetries. It is therefore believed that 
such systems may in fact be classically integrable, be- 
cause there should exist a conserved charge associated 
with each of the symmetries of the Geroch group [10]. 
Thus, these systems seem to be simple enough as to ex- 
pect that their quantization may be feasible. 

On the other hand, the existence of the Geroch sym- 
metries is on the basis of a series of solution-generating 
techniques [11] that have been developed from different 
points of view to obtain new solutions to the Einstein 
equations. Thanks to these techniques, it has been pos- 
sible to find a variety of physically interesting classical 
spacetimes with two commuting Killing fields. 

Actually, a particular family of spacetimes of this kind 
is given by the Einstein-Rosen solutions considered in 
Ref. [4] . In these solutions, the sections of constant time 
are non-compact, and the Killing fields are hypersurface 
orthogonal. In this paper, we will focus our attention on 
spacetimes which, by contrast, have closed spacelike hy- 
persurfaces and whose commuting spacelike Killing fields 
are, in general, not orthogonal. The global structure of 
the spacetimes with these properties has been studied by 
Gowdy [12], who has shown that, in this case, the sections 
of constant time must be homeomorphic to either S 1 x S 2 
(a three-handle), a three-sphere, or a three-torus (or to 



1 



a manifold covered by one of the above). Among these 
possible spatial topologies, we will limit our discussion 
exclusively to the case of a three-torus. 

A partial symmetry reduction of this Gowdy model can 
be found in Ref. [13]. Preliminary studies of its quan- 
tization, assuming the orthogonality of the two Killing 
fields, have been carried out by Berger [14]. In addition, 
Husain [15,16] has recently proposed gauge fixing con- 
ditions in the Ashtckar formulation for removing all the 
non-dynamical degrees of freedom of the model. How- 
ever, he has not performed the gauge fixing to comple- 
tion. On the other hand, Husain has not addressed the 
quantization of the system in Refs. [15,16]. Our purpose 
here is to complete the gauge fixing procedure and con- 
struct a quantum framework for the description of the 
family of Gowdy cosmologies with the spatial topology 
of a three-torus. 

The paper is organized as follows. Sec. II deals with 
the two commuting spacelike Killing field reduction of the 
Ashtekar formalism for the case of the Gowdy cosmolo- 
gies. In that section, we also present our model and dis- 
play the expressions of the first-class constraints. In Sec. 
Ill we introduce a set of gauge fixing conditions and show 
that they are well-posed and consistent. The final result 
of our gauge fixing is that we can remove all the first-class 
constraints, except for a homogeneous one. This homoge- 
neous constraint is analogous to the periodicity condition 
discovered by Gowdy [12]. The classical reduced model 
determined by our gauge fixing conditions is studied in 
Sec. IV. We prove that the phase space of the model 
can be described by using a canonical set of elementary 
variables that are all real. These are given by four func- 
tions on S 1 and two homogeneous variables. In addition, 
we explicitly obtain the metric and the reduced Hamil- 
tonian that generates the dynamical evolution. Sec. V is 
devoted to the quantization of the above reduced model 
following the canonical quantization program ellaboratcd 
by Ashtekar [2]. We first choose a representation space 
for the quantum theory and select an inner product on it 
by imposing reality conditions [2,17]. The homogeneous 
constraint of the system is then imposed a la Dirac. The 
kernel of the quantum constraint provides us with the 
Hilbert space of physical states. In Sec. VI, we deter- 
mine the form of the quantum observables of the reduced 
model and discuss the quantum evolution. Finally, Sec. 
VII contains the conclusions and some further comments. 



II. THE GOWDY MODEL 

The Gowdy universes are four-dimensional vacuum 
spacetimes with compact spacelike hypersurfaces and two 
commuting spacelike Killing fields [12]. In this paper, we 
are going to analyze only the case in which the spatial 
topology is that of a three-torus. Besides, we will restrict 
our considerations to non-degenerate Lorentzian metrics. 

Let us first introduce the Ashtekar formalism for 



Lorentzian general relativity, particularizing then to the 
Gowdy model. The Ashtekar gravitational variables can 
be taken as a densitized triad, Ef, and a SO (3) connec- 
tion, A l a , both defined on a three-manifold £ [2]. Lower 
case Latin letters from the beginning and the middle of 
the alphabet denote spatial and SO(3) indices, respec- 
tively. The SO (3) indices run from 1 to 3, and are raised 
and lowered with the metric rjij = diag(l, 1, 1). For 
Lorentzian gravity, the Poisson bracket structure is given 

by 

{A^xlE^y)} =i5 b J)8^{x-y). (2.1) 

In this formula, x and y are two generic points of S, 5 b 
is the Kronecker delta, and is the delta function on 
S. 

Provided that the metric is non-degenerate, the Ashte- 
kar variables can be expressed in terms of the triad, ef, 
and the extinsic curvature [18], k abl 

E? = e?h(e), Al=r a (e)-ik ab e b \ (2.2) 

where h = (det hat,) 1 / 2 , h ab = efe bt is the inverse three- 
metric, and T l a is the SO(3) connection compatible with 
the triad [19], 

K = -\e^ k E jb {d a El+T b ca El). (2.3) 

Here, e 1 ^ is the antisymmetric symbol, Ei a the inverse 
of the densitized triad, and T a bc the Christoffel symbols 
[18]. 

In the Ashtekar formalism, the first-class constraints 
of vacuum general relativity are [2] 

Gi = V a E? = d a E<} + cy k A{E a k = 0, (2.4) 
C a = F ab E b = 0, (2.5) 
H = e i ik Ft h E^& k =Q, (2.6) 

where F* b is the curvature of the SO (3) connection, 

F* ab = d a Ai-d b Ai + e i jk AiA*. (2.7) 

For the Gowdy universes with the topology of a three- 
torus, we can always choose spatial coordinates u, v and 9 
such that (d LJ ) a and (d„) a are the two commuting Killing 
fields. For later convenience, we will normalize the peri- 
ods of these coordinates so that 2ttuj, 2-jtv, 9 <G S 1 . All 
variables of the model must then depend only on 9 and 
the time coordinate, t. Furthermore, they have to be 
periodic in 9 £ S 1 . 

On the other hand, following Husain and Smolin [13], 
we can set equal to zero the densitized triad components 

E{ = E% = E%=E% = 0. (2.8) 

The constraints Qi, Q2, and C v are then solved by 

A\ = A 2 e = Al = Al = 0. (2.9) 
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After this symmetry reduction, and renaming A = Ag, 
E = E®, the remaining first-class constraints of the sys- 
tem can be written [15] 

G = d 9 E + J = 0, (2.10) 

C = E a L d e A L a + AJ = 0, (2.11) 
H = 2EEte L M d 6 A^ + 2AEK 

-K£ K a + K 2 = 0, (2.12) 

where a, (3 = w or v, L, M= 1 or 2, e iM is the antisym- 
metric symbol in two dimensions, and we have employed 
the notation 



± a 

j - 



AlE^ 



M ALfiP 



M 



K 
J = 



(2.13) 
(2.14) 



Thus, our Gowdy model can be described by the ten 
fields (A,E,A^,E1), which are periodic functions of 9. 
The Lorentzian symplectic structure is determined by the 
Poisson brackets 



{A(0),E(0')} = i6(0-9'), 
{A L a (0),& M (0')} = i6i6 L M 6(0 



(2.15) 
(2.16) 



6(9) being the delta function on S 1 . These fields are 
subject to the constraints (2.10-12), which will be re- 
ferred from now on as the Gauss, diffcomorphism, and 
scalar constraint, respectively. Their physical interpreta- 
tion and Poisson algebra has been discussed by Husain 
[15,16]. 

It is worth noting that the variables K£ and j£ are 
not functionally independent, for one can check that 
detK£ = det j£. Therefore, one cannot replace the 
elementary variables (A^Ef) with the eight Gauss- 
invariant quantities (K£ , j£). 

It will prove most convenient to introduce instead a 
change of phase space variables from (A^,E") to K" , 



K», K, J, and 



(2.17) 



w 



\n_q vv , (j> = arctan [ -J- I , (2.18) 
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Let us show that the above change of variables can al- 
ways be inverted in the sector of positive definite three- 
metrics. Using Eq. (2.19), relations (2.18) can be equiv- 
alently written in the form 



E^ = e w sm.^, E% = e w cos<j>. 
The definitions of x and v lead in turn to 



E^ 



vE\ + Va 7 
vEl - V^" 



v 2 E» 



v 2 E\. 



(2.21) 

(2.22) 
(2.23) 



So, given v , w € M, x > v 2 , and <fi, one can always recover 
Ef. Besides, Eqs. (2.13,14) can be seen to imply the 
identities 

K»=K-K», (2.24) 

(2.25) 



K v "= (K - 2K£)v + K»x + J^x-v\ 

that enable us to find the missing components of K£ 
from our new variables. Once K£ and E^ are known, 
we can finally obtain A^ through 



(2.26) 



E L a being the inverse of El, which can always be com- 
puted because q al3 is positive definite. 

As far as we restrict our attention to the sector of 
Lorentzian non-degenerate metrics, the variables intro- 
duced above, together with A and E, can then be re- 
garded as a set of elementary variables for our model. 
Moreover, it is easy to check from Eqs. (2.15,16) that 
they form a closed Poisson algebra. The only non- 
vanishing brackets are 

{A(9),E(0')} = {J(0),<p(0')}={K(0),w(0')} 



where 





= i6(0-0'), 


(2.27) 


{K»(0),K»(0>)} 


= -iK "(0)5(0 - 9'), 


(2.28) 


{K-(0),x(0')} = 


2ix(0)6(0 - 0'), 


(2.29) 


{K^(0),v(O')} = 


iv(0)6(0 - 0'), 


(2.30) 


{K:(0),x(0')} = 


2iv(0)5(0-0'), 


(2.31) 


{k:(0), V (o')} = 


i6(9 - 9'). 


(2.32) 



<ff> = 



EIE I3L . 



(2.19) 

From Eqs. (2.2) and (2.8), we get that q a(i = h a(3 h 2 , so 
that, for positive definite three-metrics, q a P must also be 
positive definite. Therefore, x, v and w are well-defined 
by Eqs. (2.17,18), and we must have x > 0, v,w € M, 
and x > v 2 , this last inequality coming from the fact that 



det q 



a/3 



>(x- 



> 0. 



(2.20) 



As to the variables K", K£ , K, J, and <fr, we will admit 
for the moment that they are complex. 



III. GAUGE FIXING 

We will now eliminate non-physical degrees of freedom 
from our set of phase space variables by introducing suit- 
able gauge fixing conditions. These conditions, together 
with the constraints (2.10-12), will provide us with a set 
of second-class constraints that will allow to reduce the 
model. The gauge fixing conditions that we are going to 
impose are 
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Xh = E — e = 0, 
Xg = = 0, 

Kq 



Xc 



K 



= 0, 



where 



Ko = 



K 



(3.1) 
(3.2) 

(3.3) 



(3.4) 



Here, the symbol § denotes integration over G S 1 . 

In Eq. (3.1), the time coordinate t is assumed to be 
real. This condition will be seen to fix the gauge freedom 
associated with the scalar constraint (2.12). Our gauge 
fixing is in fact equivalent to Gowdy's choice of time [12] 
(and, therefore, to that made in Ref. [16]), which can be 
expressed as E = r, t being a strictly positive time coor- 
dinate. None the less, we notice that, while in Gowdy's 
time all classical solutions present a cosmological singu- 
larity at t = [12], this singularity is driven to minus 
infinity with our choice of gauge. In this way, we allow a 
domain of definition for t that is the whole real axis. 

We will also prove that the requirement (3.2) fixes com- 
pletely the Gauss gauge of our model. From Eq. (2.18), 
this requirement implies E" = 0. Finally, we will show 
that condition (3.3) (which was employed in Ref. [16]) re- 
moves almost entirely the diffeomorphism gauge freedom. 
Some comments are in order concerning the appearance 
of Kq in this gauge fixing condition. The quantity Kq is 
known to be a classical Dirac observable of the system 
[20], in the sense that its Poisson brackets with all the 
first-class constraints (2.10-12) vanish weakly. As a con- 
sequence, Kq is a constant of motion whose value depends 
only on the particular solution that is being considered. 
This value is invariant under any gauge transformation. 
On the other hand, since K is a periodic, function of 9, it 
must admit a Fourier series of the forrro 



K= K n (t) 



And 



n— — oo 



2vr 



K n {t) 



-in& 



K 



'2tt 



(3.5) 



Condition (3.3) amounts thus to absorb all the Fourier 
coefficients K n with n ^ by means of a diffeomorphism. 

Let us now see that our gauge fixing conditions are 
well-posed. A straightforward calculation shows that 



{XH, 


j> nH} = -2m EK, 


(3.6) 


{XG: < 


f AG} = -tA, 


(3.7) 


{xc, < 


i nC} = idg(nK), 


(3.8) 



where A and n are functions on S* 1 and n is a density 
of weight —1. If n , A and dgn are different from zero, 
conditions (3.1) and (3.3) guarantee that these Poisson 
brackets never vanish for Kq ^ 0. Therefore, provided 
that Kq does not vanish, our gauge fixing conditions are 
second-class with the constraints, and hence acceptable. 

The problems found at Kq — can be obviated in the 
following sense. Using Eqs. (2.2,3), (2.8,9), and (2.13), 
it is possible to show that the variable K can be equiva- 
lently expressed in our model as 



K 



-i h k a fjh 



a/3 



(3.9) 



Here, h is again the square root of the determinant of 
the three-metric, and k a fs and h a @ denote, respectively, 
the (a/3)-components of the extrinsic curvature and the 
inverse three-metric (with a, [3 — lu or v). Then, K must 
be purely imaginary if the three-metric is positive defi- 
nite. Suppose now that Kq = 0. Since K is imaginary 
and periodic, it follows that it must vanish at least at 
one point 6q G S 1 on each section of constant time. But 
one can then easily check that all Poisson brackets of \h 
with the first-class constraints vanish at 8q, modulo such 
constraints and our gauge fixing condition. So, our gauge 
fixing is not admissible if Kq = 0. The same conclusion is 
reached if one adopts Gowdy's choice of time, E = r. As 
a consequence, the classical solutions with Kq — 0, that 
are not compatible with our gauge fixing, turn out not to 
be included in the family of cosmologies with the topol- 
ogy of a three-torus studied by Gowdy [12]. Since we are 
only interested in analyzing this family of solutions, we 
can disregard the case Kq — 0. Furthermore, we will see 
in Sec. IV that the geometry of these solutions can be 
considered invariant under a change of sign in Kq. Mak- 
ing use of this symmetry, we can set iKq G ]R + without 
loss of generality. In this way, the point Kq = will be 
driven to the boundary of our reduced phase space. Un- 
der quantization, the possible inclusion of that point will 
be physically irrelevant inasmuch as it will correspond to 
a set of measure zero in the phase space of the system. 

From now on, we will thus take Kq ^ 0. Equations 
(3.6-8) ensure then that our conditions (3.1-3) are suit- 
able to fix the scalar, Gauss and diffeomorphism gauge 
degrees of freedom. On the other hand, employing our 
gauge fixing conditions, we can solve the scalar and Gauss 
constraints to obtain the expressions for A and J as func- 
tions of the variables K£ K", x, v, and Kq: 



■.{K^dex - 2K:d e v) + e- t {K»v - K£) 



2Kq\/ x — v : 



Kq 



[{K--KZvY + {K:Y{x-v% (3.10) 



1 We assume that all the classical elementary variables, and 
in particular K, are smooth functions of 9. In fact, it suffices 
that K G C 1 ^ 1 ) for its Fourier series to converge to K at all 
points 6 G S 1 . 



J = 0. (3.11) 

This and Eqs. (3.1,2) remove the two canonically con- 
jugate pairs (A, E) and (J, <j)) as dynamical degrees of 
freedom. 
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In addition, the diffcomorphism constraint (2.11) can 
now be rewritten 



n' 



k 



dgw — 0, 



(3.12) 



where 

n' 



2{x- v 2 ) 



[K»{vdgX - 2 X dgv) - K^dg(x - V 2 )}. 

(3.13) 



Since our fields have a periodic dependence on the an- 
gular coordinate 9, IT and w can be expanded as Fourier 
series similar to that displayed for K in Eq. (3.5). For- 
mula (3.12) fixes then all the Fourier coefficients w n with 
n 7^ in terms of Kq and the Fourier coefficients of II', 



27rn: 



n ± 0. 



(3.14) 



inKg 

The coefficient wq is however left undetermined. Besides, 
integration over S 1 of Eq. (3.12) leads to the global con- 
straint 

f n' 

K = & = 0. (3.15) 



/2tt 

We recall at this point that our gauge fixing condition 
(3.3) amounts to set all the Fourier coefficients K n of 
AT, except Kq 7 equal to zero. On the other hand, Eqs. 
(2.27-32) imply that K n and w n commute under Pois- 
son brackets with the rest of our phase space variables, 
whereas 



{K n ,w m } = i6™ 



(3.16) 



We conclude in this way that our gauge fixing condition, 
together with the diffeomorphism constraint, allow us to 
eliminate the canonically conjugate pairs {K n , W- n ) with 
n^Oas physical degrees of freedom, while the homoge- 
neous components of K and w (i.e. the Fourier coeffi- 
cients Kq and wq) remain as dynamical variables. In this 
reduction process, the diffeomorphism gauge freedom is 
not totally removed, because we are still left with the ho- 
mogenous part of the diffeomorphism constraint, n o = 0. 

In order to prove that our gauge fixing procedure is 
consistent, we still have to show that the conditions (3.1- 
3) are compatible with the dynamical evolution of the 
model. This evolution is generated by the total Hamil- 
tonian constraint [2] 



H 1 



N 



-H - iN e {C — AG) — iAG 



(3.17) 



where H, C, and G are the first-class constraints (2.10- 
12), N is the densitized lapse function, N e is the only 
non-vanishing component of the shift vectorp and A is 



2 The uj and v components of the shift vector can be made 
equal to zero after the symmetry reduction (2.8,9). 



a Lagrange multiplier. Besides, N and N e are real if 
the metric is Lorentzian, and N must be different from 
zero. What we have to check then is that there exists a 
choice of densitized lapse, shift, and A such that the total 
time derivative of each of our gauge conditions vanishes. 
This total time derivative (that will be denoted by a dot) 
is given by the sum of the Poisson bracket with H T , 
{ . 7 H T }, and the partial derivative with respect to the 
explicit dependence on the time coordinate, cV After a 
careful calculation, we get that, modulo constraints and 
gauge fixing conditions, 



XH 



iN 



Kq 
v/2^ 



1 



XG = 



Xc 



= -A - iN I \Jx - v 2 K" 



27^ 



n' 



He t deN , 
K ° d B N e 



2ir 



The requirement that xh vanishes implies 



N = 



2?r 



iK 



(3.18) 

(3.19) 
(3.20) 

(3.21) 



This and xg = determine a unique A through Eq. 
(3.19). Finally, by demanding that xc = 0, we conclude 
that N e can be any function of t, N e = N e (t). We 
thus see again that the diffeomorphism gauge freedom 
has not been entirely removed, since the shift function 
N is not completely fixed. Any diffeomorphism with 
infinitesimal parameter N 6 (t) is still allowed. Note that 
such diffeomorphisms are precisely those generated by 
the only remaining constraint IIq = 0. 

For Lorentzian metrics, the shift N (t) must be real. 
On the other hand, we have commented above that K, 
and hence Kq [from Eq. (3.4)], is purely imaginary in the 
Lorentzian case. Therefore, the densitized lapse (3.21) is 
actually real. Moreover, it does not vanish for any finite 
value of Kq. This concludes the proof of consistency of 
our gauge fixing conditions with the Lorentzian evolu- 
tion. 

The final result of our gauge fixing procedure is the 
elimination of the non-dynamical fields (A, E, J, <fi) and 
Fourier coefficients (K n ,w n ) with n ^ 0. The phase 
space of the reduced model can be described by the four 
(periodic) fields K", K£, x, and v, and the two homo- 
geneous variables Kq and wq. Since all these variables 
commuted with the non-physical degrees of freedom that 
have been suppressed, the reduction of the system does 
not alter their Poisson brackets (i.e., their Poisson and 
Dirac brackets coincide). These brackets are given by 
Eqs. (2.28-32) and 



{Kq, Wq} = 



(3.22) 



which follows from Eq. (3.16). Finally, the reduced 
model is still subject to one homogeneous constraint, 
namely, H' Q — 0. 
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IV. THE REDUCED CLASSICAL MODEL 

In this section, we will analyze the dynamical evolution 
of our reduced model, discuss the reality conditions on 
the phase space variables, and obtain the expression of 
the classical metric. 



A. Classical evolution 

We have seen that the shift function N 9 can be any 
function of time, N e (t). For practical purposes, none 
the less, this shift function can always be absorbed by 
replacing the angular coordinate 9 with 



K, v = - 



9' = 0+ [ dt'N e (t'), 

J t 



(4.1) 



where to is any given time. From the periodicity of 9, it 
follows that 9' is defined as well on S 1 . Notice also that, 
for fixed t, dg> = dg. 

Since the vector fields (dgiY and (d t ) M (n = t, a) arc 
orthogonal, the dynamical evolution (with 9' kept con- 
stant) is generated in our model by the reduced Hamil- 
tonian density, Hr, which is provided by the negative 
of the momentum canonically conjugate to the variable 
chosen as time. Recalling condition (3.1), which implies 
that t = \nE, and taking into account the Lorentzian 
Poisson brackets (2.27-32), we conclude that the momen- 
tum canonically conjugate to our time variable is given 
by iAE. In this way, we arrive at Hr = —iAe*, with 
A given by Eq. (3.10). Hence, for constant 9', the time 
derivative of any of our reduced phase space variables, f, 
is 



/ = dtf + {/, & H R }. 



(4.2) 



In the following, we will assume that the shift N e (t) has 
already been absorbed in 9' and suppress the prime from 
this angular coordinate. 



B. Real phase space variables 

Our variables (x,v, K^f 7 K^, K 07 w n ) present the prob- 
lem of possessing domains of definition that are rather 
complicated in the sector of non-degenerate Lorentzian 
metrics. On the one hand, we know that x must be 
greater than v 2 , v being real. On the other hand, partic- 
ularizing to our gauge-fixed model the definition of K£ 
in Eq. (2.13) and expressions (2.2,3) for the Ashtekar 
connection, it is possible to show that 



k: 



(2xdev — vdgx) 



4(x — v 2 ) 2 
-ie* e w (x - v 2 )ik[ 



^de(x-v 2 ) 

A(x — v z )2 

-ie^e w (x- v 2 )h». 



(4.4) 



In these formulas, k„ is the extrinsic curvature, and 



w = w n e m /v27r, with w n determined by Eq. (3.14) 
for all n / 0. So, KJf and KJ are not only complex, 
but their real parts are besides functionally dependent 
on x and v in an explicitly time dependent way. We fi- 
nally recall that, for positive definite three-metrics, the 
functions w and K have to be real and purely imaginary, 
respectively. In addition, we have assumed that Kq ^ 0. 
Therefore, we must have wo € M and iKq € M + U 2R~. 

To overcome the above problems, we will now perform 
a change to a different set of elementary variables whose 
elements are all real. As a plus, these new variables will 
form a remarkably simple algebra under Poisson brackets. 

Our first step will consist in replacing x with a new 
metric variable, y, whose domain of definition is the en- 
tire real axis, 



y = In (x — v 2 ). 



(4.5) 



The condition x > v 2 ensures that y is real and well- 
defined. Using then Eqs. (2.28-32), it is not difficult to 
show that the variables 



P v = iK» (4.6) 

are the momenta canonically conjugate to y and v, 

{y(9),P v (9')} = {v(6), P v {9')} = 5(9 9'). (4.7) 

The inverse of relations (4.5,6) is 

x = e v + v 2 , 
K£ = -i{2P y + vP v ), K» - -iP v . 

On the other hand, it is obvious that 

fc = iK 



(4.8) 
(4.9) 



(4.10) 



is real and, given Eq. (3.22), canonically conjugate to 
wo, 



{w ,k } = 1. 



(4.11) 



Besides, wo and fc commute under Poisson brackets with 
y, v, P y , and P v . We have hence attained a canonical set 
of elementary variables for our reduced model. 

However, the variables P y and P v are still complex. Ac- 
tually, it follows from Eqs. (4.3,4) and (4.6) that the real 
parts of P y and P v run over the whole real axis, whereas 
their imaginary parts, T(P y ) and I(P V ), are restricted to 
be 



(4.3) 



l{P y ) = -e-i 



I(Pv) 



-e-idgy. (4.12) 
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Nevertheless, it is now easy to arrive at a set of real 
elementary variables. This can be achieved by means of 
the canonical tranformation generated by the functional 



T = w k' Q + f[(y- 2t) Pu + v Pv ] + iF, 



(4.13) 



where 



F 



*d e v (4.14) 



and p u (9), Pv{9), and k' are the new momenta. It is 
straightforward to check that kg = k' and that our con- 
figuration variables v(9) and Wq are not affected by the 
transformation. The field y(9) is however replaced with 
a new configuration variable, u(9), whose domain of def- 
inition is given again by the real axis: 



u = y — 2t. 



(4.15) 



One can see that this change of variable partly simpli- 
fies the explicit time dependence of the metric. More 
importantly, since F satisfies 



SF 
8y 



I(P, 



6F_ 

Sv 



(4.16) 



the new momenta p u and p v turn out to coincide with 
the real parts of P y and P v , respectively, 



P y = Pu + i1{Py 



P v = Pv +il(P v ). (4.17) 



Therefore, (u,p u ,v,p v ) (which are fields on S 1 ) and 
(wo,ko) provide a canonical set of real elementary vari- 
ables for our reduced model. All of these variables run 
over the whole real axis, except ko, which has to be 
non- vanishing. Besides, the model possesses one homo- 
geneous constraint, that is given by Eq. (3.15). After 
some calculus, we can express it as 



n„ 



n 



n 



o. 



-ill' = d g up u + dgvp v 



(4.18) 
(4.19) 



This constraint reproduces the periodicity condition 
found by Gowdy [12]. In fact, it was obtained from Eq. 
(3.12) by assuming that the function w is periodic. As 
remarked by Gowdy, it can also be interpreted as the 
condition that the "total field momentum be zero" [12]. 

Let us finally notice that, since the generating func- 
tional T of the above canonical transformation is explic- 
itly time dependent, the dynamical evolution of the vari- 
ables {u,p u ,v,p v ,WQ,kQ) is not generated by the Hamil- 
tonian Hr anymore. The new reduced Hamiltonian, H r , 
can be obtained from the standard formula 



Hr 



d t T. 



Using that dtT = —2 § p u + iF and Hr 
A given by Eq. (3.10), we get 



(4.20) 



-iAe*, with 



/ /, — — 



ipl + e^pl + (douf + — (devf 



(4.21) 



It is straightforward to see that this Hamiltonian is 
bounded from above (below) for positive (negative) val- 
ues of k . 



C. The metric 

We will now obtain the expression of the classical met- 
ric that results from our gauge fixing conditions for the 
Gowdy model. Using Eqs. (2.2), (2.8), (2.17-19), and 
the definition of u, we get that the only non- vanishing 
components of the three-metric h a b are 



ha 



e 2w e : 



9af3, 



(4.22) 

9uu = 1, g U u = -v, g vv = e 2t e u + v 2 . (4.23) 

The lapse function, N — N (det Ziab) 1 ^ 2 , can then be 
found from Eq. (3.21), 



N = 



'2-k 



e t e w e- 



(4.24) 



Introducing now the change of time coordinate 

T=^L*, Tel+, (4.25) 



and remembering that the shift function N 9 has been 
absorbed in the angular coordinate 9, we arrive at a line 
element of the form 

ds 2 = e 2w e% (-dT 2 + d9 2 ) 

+ (gujuidiu 2 + 2g UJl/ dujdi' + g vv dv 2 \ (4.26) 

where, in terms of the positive time coordinate T, g vl/ 
reads 



9 ,. = grV 



(4.27) 

On the other hand, Eq. (3.14) determines w to be 

oo „ 

w ^ IE. 



w 



y 

inkn 

= -oo,^0 U 



n e in6 



(4.28) 



with n„ the Fourier coefficients for II, and Uq — be- 
cause of constraint (4.18). We have thus succeded in 
writing the metric of our model in terms of the dynami- 
cal variables u, p u , v, p v , wq and fco- 

It is easy to check that the classical geometries de- 
scribed by Eqs. (4.26-28) are in fact invariant under a 
flip of sign in the momenta p u , p v , and kg. a We can take 



This change of sign in the momenta, while keeping unal- 
tered the configuration variables, can be interpreted as a time 
reversal. 
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advantage of this symmetry to fix, e.g., the sign of fco, 
without eliminating from our considerations any of the 
geometries that are allowed for our model. We will hence 
restrict fco to be positive from now on: fco G St + ■ 

Given this restriction, it is convenient to replace 
( wo, fco) with a new pair of canonical variables, (bo, Co), 
whose respective domains of definition are the whole real 
axis, 



b = k w , 



c = In fc - 



Note that 



{bo, c } = 1, 



(4.29) 



(4.30) 



and that the point fco = 0, excluded from our phase space, 
has now been driven to the boundary of the domain of Co 
(namely, to c = -co). Notice also that, after performing 
the change (4.29), the Hamiltonian (4.21) turns out to be 
analytic in all the elementary fields and variables of our 
reduced model. 

To close this section, we will derive a formula equiva- 
lent to Eq. (4.28) that may be more useful in practice 
to compute the metric function w. A straightforward 
calculation leads to 



ri = {n, J> H r } = d e H r , 



w = {w , f H r } = -— f H r . 



(4.31) 
(4.32) 



It is worth remarking that the first of these equations 
implies that the constraint (4.18) is preserved by the dy- 
namical evolution of our model. From Eqs. (4.31,32), we 
also get the following relations among wo and the Fourier 
coefficients of II and if r (n„ and if™, respectively): 

— = / dt'H?(t') + d n , n^O, (4.33) 
in Jo 

wo = / dt'H?(t') + d , (4.34) 

fc(D Jo 

the dnS being constants for all integers n. Using these 
relations, it is not difficult to check that Eq. (4.28) can 
be rewritten as 



w(t, 0)-- 



fco 
+D, 



[ dt'H r (t',e = o)+ [ de'u(t,e' 

Jo Jo 

(4.35) 



where D is an undetermined real constant. 

V. QUANTIZATION 

We have seen that the four fields (u,p u ,v,p v ) and 
the two homogeneous degrees of freedom (&o, c o) form a 
set of real elementary variables for our reduced model. 



They are only restricted in that they are subject to 
the constraint (4.18). Since all our elementary fields 
are defined on S 1 , we can expand them as Fourier se- 
ries. The corresponding Fourier coefficients will be called 
(u n ,p™, v n ,p™). These coefficients, together with (bo, Co), 
provide then an infinite set of homogeneous, elementary 
variables for our system. From the fact that (p u ,p v ) are 
the momenta canonically conjugate to (u, v), we arrive 
at the following Poisson bracket structure: 

{u n ,p™} = {v n ,p™} = 5™ n , {6o,c } = l, (5.1) 

the rest of brackets being equal to zero. On the other 
hand, the reality conditions on our fields (u,p u , v,p v ) and 
variables (bo,Co) imply that 



bo = bo, co = c , 



(5.2) 



where g = u or v and the bar denotes complex conjuga- 
tion. 

In order to quantize the system, we will follow the 
canonical quantization program put forward by Ashtckar 
[2]. We will first represent our elementary variables by 
linear operators acting on an auxiliary vector space. An 
inner product will be selected on this space by imposing 
the reality conditions (5.2) as adjointness relations be- 
tween our operators [2,17]. We will then represent the 
constraint (4.18) as an operator, IIo, and impose it on 
our quantum theory. From the kernel of IIo and the in- 
ner product introduced on the auxiliary representation 
space, we will construct the Hilbert space of physical 
states, Hp. Finally, we will identify the quantum observ- 
ables of our reduced model. By quantum observable, we 
mean any operator that has a well-defined action on TL p 
and, therefore, commutes with the only constraint of the 
model, IIo- In the rest of this section, we will implement 
this quantization program, except for the discussion of 
the quantum observables, that will be postponed to Sec. 
VI. At some stages of our analysis, we will proceed only 
in a formal way; thus, our main concern is to show how 
the quantization process can be carried out, rather than 
obtain explicitly a rigorous and complete quantum the- 
ory for our model. 

We begin by choosing as our auxiliary representation 
space the vector space of analytic functionals \f of the 
set of variables Q, = (co,u n ,v n ) [n — 0, ±1, ...]. On this 
space, we can represent our elementary variables by the 
operators 



dco' 



(5.3) 



5n* = 9n *, P g n * = -*-5— , (5-4) 

where g stands again for iioru (and we have set h = 1). 
Notice that the commutators of these operators repro- 
duce the classical Poisson algebra (5.1) with the due fac- 
tor of i: 
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[bo, c ] 



[9n,P~„ 



(5.5) 



The operators of our quantum theory are then given by 
(possibly infinite) sums of products of the elementary 
operators (5.3,4) and c-number operators [2]. 

Using the reality conditions (5.2), we can select an in- 
ner product on our auxiliary representation space in the 
following way. We first adopt the ansatz 



< $ 



<■>>=/ 



dnAdQp{Q,n) tf(fi), (5.6) 



with 



dfl A dfl= -dc A dco x 

[[ ( l -du n A du n ) JJ ( ^du m A df 



(5.7) 



and p a certain positive integration measure. The re- 
quirement that the reality conditions (5.2) be realized as 
adjointncss relations determines then the measure p up 
to an overall positive constant, 

p = S(co - Co) J| [<5(u_„ - u^) S(v- n - WO] • (5-8) 



In this formula, we have employed the notation 

5(z)8(z) = 5(x)5(y), 5{z -z) = 5(y), (5.9) 

for z = x + iy any complex variable, and x,y e M. 

We still have to impose the homogeneous constraint 
(4.18) quantum mechanically. In order to do it, let us 
represent the phase space variable n by the operator 



n = j2s(x:+x:), 

s=l 

where (for g = u or v) 



(5.10) 



(5.11) 



We note that, with the factor ordering chosen in Eq. 
(5.11), IIo is, at least formally, self-adjoint on the auxil- 
iary Hilbert space determined by the inner product (5.6- 
8).' 

The physical states of our quantum theory are those 
annihilated by the constraint operator IIo. In our auxil- 
iary representation space, on the other hand, all quantum 
states ^ can be expressed as (possibly infinite) sums of 
functionals of the form 



n 



(5.12) 



where a = (i n ,j n ) [n = 0, ±1,...], and k, i n , and j n 
are non-negative integers. From Eqs. (5.10,11), it is 
easy to check that the functionals P(k,a) are (generalized) 
eigenfunctions of n , 



floP(k,*) =N(a)P {k ^. (5.13) 
Here, N(a) is given by 

oo 



and is equal to the sum of indices of all factors {u ni v n ) 
(including multiplicities) appearing in P(k,a)- F° r an Y 
quantum state ^ — Y^, a (k,a)P(k,a), with a.(k,o) some com- 
plex constants, we therefore get 



ft * = a (k,a)N{(T)P {kta) . 



(5.15) 



The uniqueness of the power series of the zero functional 
implies then that \1/ is annihilated by IIo if and only if 
it is a (generally infinite) linear combination of function- 
als P(k,a) whose total indices N(a) vanish. The ana- 
lytic functionals that satisfy this condition form a com- 
plex vector space, whose Hilbert completion with respect 
to the inner product (5.6-8) finally provides us with the 
Hilbert space of physical states, H p . 

It is easy to see that this Hilbert space is actually in- 
finite dimensional. For the sake of an example, let us 
display an infinite set of states in H p , namely, 



exp ■ 



C 

1 

I 



n 

,s>0 



C 2 



1 



2nA,B, 



E 



U S U- 



V S V- 



Bi 



(5.16) 



where C, A s , and B s (s = 0,1,...) are real constants, 
but otherwise unrestricted, and P(k,<r') is any polynomial 
of the form (5.12) (i.e, the set a' contains only a finite 
number of non- vanishing elements in this case) such that 
N(a') = 0. All these states are analytic functionals of 
il, belong to the kernel of ITo and can be checked to be 
normalizable with respect to the inner product (5.6-8). 



VI. QUANTUM OBSERVABLES AND 
HAMILTONIAN 



We turn now to the task of finding the quantum ob- 
servables of our model. These are the operators with a 
well-defined action on the Hilbert space of physical states, 
Hp. As they leave TC P invariant, they must (weakly) com- 
mute with the quantum constraint of the model, Ho. 

In our quantum theory, on the other hand, all opera- 
tors are supposed to be constructed from the elementary 
operators (5.3,4) by taking (suitable limits of) sums and 
products. In particular, it should be possible to obtain all 
quantum obscrvables from (infinite) linear combinations 
of operators of the form 
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*W) = # b k 2 II [<" vt (Pu n ) 1 - (AT")™-] , (6.1) 

n 

k = (ki,k 2 ) and T = (i n ,j n , l n , m n ) [n = 0, ±1, ...] being 
two sets of non-negative integers. 

A straightforward calculation shows that 



(*,r)J 



N(T)P { 



where 



N(T)=^Ts{i s +j s +l s 



fc,r). 



+ m a - i- a 



(6.2) 



s=l 



-i- s - l-s - m- 



(6.3) 



is [similarly to N(a) in Eq. (5.14) ] the sum of indices 
of all the elementary operators (counting multiplicities) 



that form Pi 



k,T)- 



Therefore, P(k,r) commutes with Iln 



if 



and only if its total index iV(r) vanishes. Furthermore, 
from Eq. (6.2) and our comments above, it is possible 
to show that the quantum observables of our theory can 
always be expressed!!! as linear combinations (including 
the limit of infinite sums) of operators P^ p) verifying 
N(T) = 0. 

A possible way to attain observables is the following. 
We first define 



m= E 



pm= E 



'2it 



(6.4) 



Here, g = u or v, and 9 S S 1 must be regarded as a 
parameter. Suppose next that g(co,&o) is an analytic 
function of Co and 6o, and f[u(6),v(6)] a functional of 
only u{9) and v(9) (and perhaps of their derivatives with 
respect to 9) which is analytic in these fields. One can 
then check that, for all non-negative integers I and m, 
the operators 



g{c ,b ) f f[u,v](pu) 1 (p v 



(6.5) 



can be written as (infinite) sums of operators P(k,r) with 
N(T) = 0, so that they provide, in general, quantum 
observables for our model. 

So far, we have not discussed the dynamics of our quan- 
tum system. In order to do it, we must first allow the 
physical states and quantum observables to depend on a 
dynamical parameter t S M. The quantum evolution is 
then dictated by the Schrodingcr equation 



(6.6) 



4 Up to additive terms of the form XTlo (X being a generic 
operator), which vanish modulo the constraint IIo. 



where Hj is, by assumption, a self-adjoint observable 
that represents the classical reduced Hamiltonian inte- 
grated over S 1 , that is, the generator of the classical evo- 
lution, § H r . 

The self-adjointness of Hj implies that the quantum 
evolution is unitary, i.e., that it preserves the norm of 
the physical states of the theory. This is equivalent to 
say that the integration of the Schrodingcr equation (6.6) 
leads then to a unitary evolution operator, U(t, 0), such 
that = U(t, 0)*, * being the initial physical state. 
In terms of the quantum Hamiltonian, the evolution op- 
erator adopts the expression [21] 



U{t,0) =V 



exp 



dt'H?(t'] 



(6.7) 



where V stands for the time ordering 



V[Hj(h) . . . Hj{t n )\ = E #r r (Vi)) • • ■ Hj(t v{n) ) x 
(v) 

@(*tj(1) _ V2)) • ■ ■ ©(*»?(n-l) _ *J7(n))- (6-8) 

Here, is the Heaviside function and r\ any permutation 
of the indices 1,. . ., n. 

From the Schrodingcr equation and the self-adjointness 
of the quantum Hamiltonian, we also arrive at the fol- 
lowing evolution for the matrix elements of any quantum 
observable O(t): 



i- <$(*), 0(t)*(t) >= 



< $ 



(t), ([6{t),Hj{t)]+id t 6{t)) *(t) >, (6.9) 



with d t O{t) the derivative of 0(t) with respect to its 
explicit dependence on the parameter t. We will then say 
that an observable 0(t) represents a constant of motion 
if it satisfies 



[d(t),H^(t)] + td t d(t) = 0, 



(6.10) 



so that all its matrix elements are constant in the quan- 
tum evolution. In this sense, it is worth pointing out 
that, given any quantum observable, O, that is explicitly 
^-independent, one can generally obtain another observ- 
able that represents a constant of motion, namely, 



6'{t) = U(t,0)OU' 1 {t,0), 



(6.11) 



where C/ -1 (£,0) is the inverse of the evolution operator. 

We have thus seen that, in order to arrive at a unitary 
quantum evolution and essentially complete our quanti- 
zation, we are only left with the problem of finding a self- 
adjoint observable to represent the (integrated) classical 
reduced Hamiltonian of the model. A quantum Hamil- 
tonian that, at least formally, satisfies these conditions 



H „ 



-X 



„2( 



Y, 



.12) 
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X 



'2ne 



A(p u f + -e^(dev) 2 



E 



^p; n - E T e 



-n—m v nVm 



m— — oo 



(6.13) 



27re- C0 

OO 

E 



16 



{d e uf + e*{ Pv f 



16 



n 2 u n u_ n + e ~n- m PvPv 



In the above formulas, 



-ind 



2ir 



(6.14) 



(6.15) 



and u, p„, and p„ are the operators defined in Eq. (6.4). 
It is clear that this quantum Hamiltonian commutes with 
H , because it is a linear combination of operators of 
the form (6.5). That this Hamiltonian is formally self- 
adjoint follows from the fact that it is given by a sum 
of products of commuting operators, as well as from the 
reality conditions = c , g\ — g- n and (]3™)^ = pg' n 
(g = u or v). Finally, notice that Hj inherits an explicit 
dependence on the parameter t from the time dependence 
of the classical Hamiltonian (4.21). 

To prove that the Hamiltonian (6.12-14) is in fact a 
self-adjoint observable, it would actually suffice to show 
that it is densely defined on the Hilbert space of phys- 
ical states, Hp. From our discussion above, this would 
guarantee that Hj is a symmetric observable. That this 
Hamiltonian is self-adjoint (or, strictly speaking, that it 
admits a self-adjoint extension) would then be a conse- 
quence of the fact that there exists a conjugation C on H p 
which leaves the domain of Hj invariant and commutes 
with it [22]. We remind that a conjugation C : H p —* H p 
is an anti-linear, norm-preserving map whose square is 
the identity. It is not difficult to check that a map on TL p 
that satisfies the properties of a conjugation and com- 
mutes with our quantum Hamiltonian is 



C*(fi) = *(C(fi)), 



(6.16) 



where is any physical state and the action of C on the 
set of elementary variables O is given by 



C c = c , C u r , 



Cv n 



(6.17) 



with n = 0, ±1, ... 

It could also happen that, instead on Hp, the Hamil- 
tonian (6.12-14) admitted a self-adjoint extension only 
on a sufficiently large Hilbert subspace H p C H p . In 
that case, one could still try to restrict all considera- 
tions to that subspace in a consistent way, and regard 



H p as the true Hilbert space of physical states. Other- 
wise, one would have to replace the operator (6.12-14) 
with a different quantum Hamiltonian that turned out 
to be physically acceptable in our model. Finally, if no 
such Hamiltonian could be found (and one insisted in ar- 
riving at a unitary quantum evolution), one would have 
to start the quantization over again, changing any of the 
choices that are available in the construction of the quan- 
tum theory, like, e.g., the set of elementary operators or 
their representation. 



VII. CONCLUSIONS AND FURTHER 
COMMENTS 

Starting with the Ashtekar formalism for Lorentzian 
general relativity in vacuum and restricting our atten- 
tion to the sector of non-degenerate metrics, we have 
discussed the structure of the reduced phase space and 
the quantization of the family of Gowdy universes whose 
spatial topology is that of a three-torus. 

We have first removed non-physical degrees of freedom 
by means of a gauge fixing procedure. The gauge fix- 
ing conditions imposed, together with the first-class con- 
straints of the model, have been shown to form a set of 
second-class constraints that allow the reduction of the 
system. In this way, we have been able to eliminate all 
the constraints of the model except for one homogeneous 
constraint, Ho = 0. This constraint is the analogue of the 
periodicity condition studied by Gowdy [12], and gen- 
erates the diffeomorphisms, with spatially constant in- 
finitesimal parameters, of the angular coordinate 9 that 
does not correspond to a Killing field of the spacetime. 

The choice of time that we have adopted is equivalent 
to that employed by Gowdy [12,16]. We have got rid 
of the Gauss constraints and of the diffeomorphism con- 
straints of the coordinates associated with Killing fields 
by requiring that some components of the densitized triad 
vanish [see Eqs. (2.8), (2.18), and (3.2)]. Finally, the In- 
coordinate diffeomorphism gauge freedom has been used 
to set the variable K [given by Eq. (3.9)] equal to its 
mean value on each surface of constant time [16], i.e., to 
KqI \p2/K. This quantity is known to be a constant of mo- 
tion of the model. We have then shown that the classical 
geometries with Kq = are not included in the family 
of cosmological solutions considered by Gowdy. Besides, 
provided that K is different from zero, our gauge fixing 
conditions are consistent and well-posed. 

We have found a canonical set of real elementary vari- 
ables for the phase space of our reduced model. This 
set is formed by the four fields u(6), p u (&), v(9), and 
p v (0), and by the two homogeneous variables wq and ko- 
The reduced model is still subject to the homogeneous 
constraint Ho = 0. On the other hand, the exclusion of 
the solutions with A'o = implies that fco € M + U M~ . 
Making use of the fact that the classical geometries are 
invariant under a change of sign in the momenta p u (0), 
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p v (0), and fco, we have none the less been able to restrict 
all considerations to the case fco € M + without loss of 
generality. In order to attain a canonical set of elemen- 
tary variables whose respective domains of definition are 
the entire real axis, we have then replaced wo and fco with 
a new canonical pair of variables, (60, cq). 

In addition, we have obtained the explicit expression 
for the classical metric of the Gowdy spacetimes and de- 
termined the reduced Hamiltonian, H r , that generates 
the dynamical evolution in our gauge- fixed model. This 
Hamiltonian presents an explicit dependence on the time 
coordinate, so that the reduced system is not conserva- 
tive. 

Since the fields u(9), p u (6), v(0), and p v (9) are peri- 
odic functions of 9, we can expand them as Fourier series. 
The Fourier coefficients {u n ,p~ n ) and (v n ,p~ n ) turn out 
to be canonically conjugate pairs of homogeneous vari- 
ables. Employing these Fourier coefficients and (bo,co) 
as elementary variables, we have proceeded to quantize 
our model following the canonical program ellaborated 
by Ashtekar [2]. We have first represented the variables 
60, Co, u n , p™, v n , and (n = 0, ±1, ...) as elementary 
linear operators acting on the vector space of analytic 
functionals of Co, u n , and v n . A unique inner product 
has been selected on this space by demanding that the 
complex conjugation relations (5.2) (our reality condi- 
tions) are realized quantum mechanically as adjointness 
relations. We have then represented the homogeneous 
constraint of our reduced model by a linear operator, ITo, 
and determined the quantum states that are annihilated 
by it. These states, together with the inner product se- 
lected by the reality conditions, have provided us with 
the Hilbert space of physical states, 7i p . 

The quantum observables of the reduced model are the 
operators that have a well-defined action on H p . In our 
quantum theory, on the other hand, a generic observable 
should always be given by a suitable (possibly infinite) 
sum of products of elementary operators. Using this fact, 
we have been able to obtain the general form of the quan- 
tum observables. 

We have finally introduced a dynamical evolution in 
our system by imposing a Schrodinger equation with 
quantum Hamiltonian, Hj ', representing the classical re- 
duced Hamiltonian integrated over S 1 . If one requires 
that the quantum evolution be unitary, the Hamiltonian 
Hj must be a self-adjoint observable. We have found an 
operator Hj that, at least formally, satisfies these condi- 
tions. Also discussed are other still available, alternative 
possibilities to obtain a Hamiltonian which would really 
be well-defined and self-adjoint. 

In analyzing the structure of the phase space of our 
reduced model, wc have restricted the variable fco to be 
positive by taking advantage of the symmetry of the clas- 
sical geometries under a change of sign in the momenta 
p u (0), Pv(9)j and fco, a transformation that can be re- 
garded as a time reversal. Had we not imposed this re- 
striction, we should have split the phase space into two 



disconnected parts: one for fco > 0, and the other for 
fco < 0. Replacing definiton (4.29) with 

b = k w , c = In (-fc ) (7.1) 

in the sector of negative values of fco and repeating our 
quantization procedure, we would have then arrived at a 
quantum theory whose physical Hilbert space would be 
given by the direct sum of two copies of the Hilbert space 
of physical states constructed for fc € M + . Nevertheless, 
any of these two copies would actually provide us with 
an irreducible representation of the model as far as we 
do not allow time reversal operations. 

On the other hand, although we have considered the 
invariance of the Gowdy geometries under the transfor- 
mations generated by the diffeomorphism constraints, we 
have in fact not discussed the possible symmetries under 
global diffeomorphisms that cannot be connected with 
the identity transformation. Diffeomorphisms of this 
kind which are compatible with the form of the met- 
ric (4.26) arc given, e.g., by a change of orientation in 
one of the angular coordinates, or by an interchange of 
the coordinates that correspond to Killing fields of the 
spacetime, i.e. u> and v. In this sense, the point of view 
that we have adopted is that the coordinates ui and v 
are physically distinguishable and the orientation of all 
spatial coordinates fixed once and for all. 

During the completion of this paper, we have become 
aware of an independent work by A. Ashtekar and M. 
Pierri [23] , who also study the quantization of the family 
of Gowdy universes with the spatial topology of a three- 
torus. In that work, the discussion has none the less been 
restricted to the case in which the two commuting Killing 
fields of the model are hypersurface orthogonal. 
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